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It is an interesting question whether classical gamma ray bursts repeat, espe¬ 
cially because of the potential implication on the distance of the sources. 
Gamma ray bursts were first discovered [1] by Vela satellites in 1967 the 
same year the first radio pulsar was discovered. Pulsars are understood to 
be rotating neutron stars in essence, but the origin of gamma-ray bursts 
is a complete mystery. If they are from the (extended) Galactic halo, the 
isotropic emission energy will be ~ lO^^ergs. If they are from cosmological 
distances, the isotropic emission energy will be ~ lO^^ergs. For example, the 
bursts from cosmological distances can not be repetitive because the enor¬ 
mous energy output would most likely be accompanied by the consumption 
of the sources. 

Gurrently, we do not have dynamic clues toward a more meaningful iden¬ 
tification or rejection of repetitions, and the controversy lies in whether there 
is statistical evidence of repeating classical gamma ray bursts in the BATSE 
catalog. In order to address the problem, we ask ourselves the following 
questions. 

1. Can we define a quantity which measures the probability for a given pair 
of GRB events to have come from the same position of the sky? 

2. Can we find a statistic which can gauge the significance of matched pairs? 

3. Is there a statistically significant excess of matched pairs in the BATSE 
catalog? 

As an answer, we construct the ‘match probability’ of a given pair of events 
based on measurement error distribution functions, build statistics using 
‘the match probability’, and analyse the publically available BATSE IB and 
2B data sets. 

Our statistics is built based on the simple observation that the repetitive 
bursts from the same source form an error distribution in the position space. 
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An obvious but important aspect of it is that it is highly localised (around 
the true source position). Gamma ray bursts in the BATSE catalog are 
known to be largely isotropic, and the BATSE catalog can be considered as 
a mixture of error distributions and a background isotropic distribution. An 
isotropic distribution is characterized by the complete lack of localization. 
The drastic difference in localization behavior of the two types of distribu¬ 
tions is the underlying basis of our burst pair match probability statistics. 

Given a pair of bursts, the probability that they have come from the same 
source depends on the joint probability density of the error functions of the 
events /(xi) and /(X2). 


P 12 = [ /l(xi)/2(x2) = [ /l(x - Xic)/2(x - X2c)d2x , 

Js^ Js^ 


( 1 ) 


where xic and X 2 c are the center positions. If we assume gaussian error 
distributions (in 3?^, which is good enough here), 


P 12 = 


2710-^2 




( 2 ) 


where (T ^2 = + ^ 2 - define the ‘match probability’ of the pair by 

integrating pi 2 from the measured distance ri 2 = |xi — X 2 I to infinity. 


P12 = P(ri2,ui2) = -\r r-dr = e-"?2/2-f2 (3) 

Now, we look at the distribution of pairs in P-space. For an error distri¬ 
bution in position space, the pair distribution in P-space is homogeneous. 
In other words, the number of pairs is a constant function of P. Since 
P G [0,1], the constant is the total number of pairs. By normalizing by the 
total number of pairs, we define the pair density distribution g{P). 


5err(E) = l. ; E G [0, 1] 


(4) 


For the background isotropic distribution, 

2 

5iso(-P) ot ^ • (5) 

(The number of pairs at a distance r oc rdr and r and P are related by Equa¬ 
tion (|3|.) It is not surprising that giso{P) is highly concentrated toward small 
P because isotropic distribution has many more pairs at large separations 
in comparison to an error distribution of the same number of bursts. 

u << I for most of the bursts in the BATSE catalog, and ( 7 err(E) > 
gisoiP) unless P is very small. Thus, we define the following statistical quan¬ 
tities. 


^Qa — Qaip) ^ QaijP) Piso j — 0, 1, ... 


( 6 ) 
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i<j 


( 7 ) 


where 0(x) = 1 if x > 0 and 0(x) = 0 otherwise. The signal should be 
found in A.Qa, and the background isotropic distributions for < Qa{p) >iso 
are generated by randomizing the burst positions while keeping the error 
functions of the data. In the limit of a large sample, we can calculate AQa.{p). 

= T. ( 8 ) 

matched pairs 

Therefore, we can calculate the expected signal to noise ratio (S/N) by com¬ 
paring AQat\i{p) to the RMS fluctuations about < Qa{p) >iso- We find that 
gi(0) has the highest S/N. For the fraction of the matched pairs (of bursts 
G error distribution) of 10“^, the S/N of gi(0) 3.0 for 485 bursts (exclu¬ 

sive of the events without error estimations) of BATSE 2B. We emphasize 
that our ‘match probability’ test of the BATSE catalogs for repetition is at 
the match level of ~ 10“^. For example, if five bursts in the BATSE IB are 
from a repeater as suggested by Wang and Lingenfelder [6], the signal can 
not be seen statistically because the required sensitivity is 3 x 10“^. The 
analyses of gi(0) statistics of the BATSE data are shown in Table 1. 


TABLE I 

The number of matches and anti-matches measured with the Qi(0) statistic 
with 1-a error bars determined from the RMS deviation from the mean of 
10,000 simulated BATSE catalogs. The significance refers to the fraction of 
the 10,000 simulated catalogs with larger values of Qi(0) than the real data. 


Sample: 

All 

fTstat < 9° 

O-stat < 4° 

Exposure Cut 

-^match IB 

13 ±25 

27± 14 

1±7 

3± 12 

-^antipode IB 

54± 25 

24± 14 

17 ±7 

9± 12 

match signif. IB 

29.3 % 

2.65% 

43.36% 

36.96% 

antipode signif. IB 

2.42% 

4.54% 

1.32% 

21.92% 

T/match 2B 

-23 ± 45 

0±26 

-13± 13 

-9 ±20 

A^antipode 2B 

31 ±45 

39 ±26 

29± 13 

1 ± 20 

match signif. 2B 

68.31% 

48.33% 

82.39% 

65.61% 

antipode signif. 2B 

23.34% 

7.84% 

1.90% 

46.72% 


We conclude that there is no signal of recurrent bursts in the BATSE 
data at the match level of 10“^. For more details, see [2]. 
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Now, where does our statistic lead to in regard to the current controversy 
on the repetition of classical gamma ray bursts? The outcome of devising a 
new statistic may not be of much value if it merely provides a third statis¬ 
tic and another opinion on the issue. It turns out that we can understand 
both two-point correlation [4] and nearest-neighbor statistics [3] from ‘match 
probability’ statistics because of the fuller position measurement informa¬ 
tion contained in it and moreover we can estimate the S/N of the statistics. 
Two point correlation function is directly related to Qo{p) and the S/N at 
the optimal pop is somewhat worse than the true optimal case of (5i(0). One 
should keep in mind that the optimal value pop can be found only by embed¬ 
ding the two point correlation function in the match probability statistics. 
In order to understand the nearest neighbour statistics (NNS), we studied 
the cases of 60 matched pairs. We find that NNS has slightly better sen¬ 
sitivity than Qi(0) when each source bursted only twice. For multiplicities 
bigger than two, the S/N of NNS degrades dramatically and is far worse 
than that of the two point correlation statistics. Therefore, the conjecture 
by Quashnock and Lamb [3] based on nearest neighbour correlations that 
there are recurrent bursts with high multiplicities in the BATSE catalog is 
not supported by more sensitive tests. NNS also has a peculiar disadvantage 
that the sensitivity does not improve as the size of the sample grows. 

Finally, we would like to comment on the ‘signal’ Quashnock and Lamb 
[3] and Narayan and Piran [4] have seen. Is the ‘signal’ identified in our 
analyses? Yes: In the sample Ustat < 9° ( without exposure corrections), the 
curve Qi{p) follows that of Qijso and then deviates from it at around p = 0.6 
to follow the (theoretical) curve of signal for 10“^ matches. Speaking in terms 
of g{P), there is a fluctuation (depletion and excess) in g{P) at around P = 
0.6 compared to that of giso- Is the ‘signal’ consistent with the expectation 
from repeating bursts? Clearly, it is not. What is this ‘signal’ a signal of? We 
do not know. It could be due to selection effects [4, 5] or many other reasons 
as yetuntested. But it is also consistent with statistical fluctuations. When 
we make exposure corrections using the trigger efficiencies provided in the 
BATSE catalog, the ‘signal’ disappears. We note that the events excluded 
by this cut are mostly faint bursts. 
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